NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 





^ MAR19at 

i iminfi 


By 

BruM J. Blanchard 
MarshaHJ. McFarlMd 
Sidney Thais 
John Q. Richlar 
Ramola Sanaing Cantar 
Texas AS M UnIvarsIty 
CoNaga Station, Texas 77843 


February 1961 


Prepared lor 

Qoddatd Space FNgM Center 
Qreenbalt. Maryland 80771 

Contract No.: N808193 


Rnal Rapoit R8C*3622*4 


(MASA-CU- 164142 ) COaHELATlOM Of SPACECBAFT 
PASSIVE MlCaOitAVE SVSTEfl DATA WITb SOIL 
nuiSIUHE INDICES (API) fioal Eepoct (Texas 
A 6 rt Univ.) 49 p HC AUS/rtf A 01 CSCL J 8 H 


GJ/43 


CORRELATION OF SPACECRAH’ PASSIVE MICROWAVE 
CYSTEM DATA WITH SOIL MOISTURE INDICES (API) 


HUI-214116 


Uucias 

415BJ 


Final Report 3622-4 


Correlation of Spacecraft Passive Microwave 
System Data with Soil Moisture Indices (API) 


By 


Bruce J. Blanchard 
Marshall J. McFarland 
Sidney The is 
John G. Richter 


Prepared for 

Goddard Spact Flight Center 
Greenbelt, Maryland 80721 


Contract No. : NSG-5193 


TABLE OF CONTENTS 


Page 

LIST OF FIGURES 11 

LIST OF TABLES v 

1.0 INTRODUCTION. 1 

1.1 Necessity of This Research 1 

1.2 Background 3 

1.2.1 Soil Moisture Estimation 3 

1.2.2 Plant Stress Estimation . . . 4 

1.3 Objectives. 7 

2.0 SITE DESCRIPTION 8 

2.1 Data Available 8 

2.2 Cropping Practices 13 

3.0 NIMBUS 5 ESMR DATA 20 

4.0 TECHNIQUES OF DATA HANDLING 22 

4.1 The Soil Moisture Model 22 

4.2 Emissivlty Model 22 

4.3 Stress Degree Index Model ...... 24 

5.0 ANALYSIS AND DISCUSSION "JO 

5.1 Correlation Between API and Emissivlty 30 

5.2 Free Water After Storm Frontal Passage 38 

5.3 Drought Conditions and Winter Wheat Yield 60 

6.0 CONCLUSION. 82 

REFERENCES 85 

APPENDIX A 88 


i 






LIST OF FIGURES 

Figure Page 

1 Southern great plains area used as a basis for calibration 

of ESMR 9 

2 Gridded area used by McFarland In preliminary study and as 

training site grid 10 

3 Soils map of the preliminary study area . 11 

4 Percentages of preliminary study area devoted to winter wheat 

and total croplands by county. The first number is percent 
winter wheat. The second Is percent croplands Including winter 
wheat 14 

5 A graymap representing the percentage of each grid point covered 

by cropland 15 

6 A graymap depicting the percentage of each grid point covered 

by continuously cropped dryalnd winter wheat 16 

7 A graymap representing the 1974 continuously cropped dryland 

winter wheat yield, in bushels per acre, for each grid point .. . 17 

8 Number of weeks that the phenological development of wheat lags 

that of the Texas crop 19 

9 ESMR ground resolution as a function of scan angle from nadir. . . 21 

10 A plot of the crop susceptibility curve developed for this study . 27 

11 A plot of the mean dally maximum and minimum temperatures used 

in model 2 29 

12 Fall (August 12-November 1) correlation coefficients for each 

grid point (values are multiplied by -100) 31 

13 Spring (February 28-April 15) correlation coefficients for each 

grid point (values are multiplied by -100) 33 

14 Fall scatter plot for six winter wheat grid points 34 

15 Fall scatter plot for six non-wheat croplands grid points. .... 35 

16 Spring scatter plot for six winter wheat grid points ....... 36 

17 Spring scatter plot for six non-wheat croplands grid points. ... 37 

18 API - Date: 3/15/74 39 

19 Emissivity - Date: 3/15/74. ........... . . 40 


ii 


LIST OF FIGURES (continued) 


Figure Page 

20 API - Date: 4/22/74 41 

21 Emissivlty - Date: 4/22/74 42 

22 API - Date: 8/25/74 . . . . 43 

23 Emissivlty - Date: 8/25/74 44 

24 API - Date: 9/22/74 45 

25 Emissivlty - Date: 9/22/74 46 

26 API - Date: 3/18/75 47 

27 Emissivlty - Date: 3/18/75 48 

28 Percent of each grid area covered by winter wheat (derived by 

state agricultural statistics) 51 

29 Linear correlation coefficients as an indicator of scatter in 
the data between emissivlty and API for Fall (August 12- 

November 1) 52 

30 Linear correlation coefficients as an indicator of scatter in the 
data between emissivlty and API for Winter (November 2- 

February 27) 53 

31 Linear correlation coefficients as an Indicator of scatter in the 
data between emissivlty and API for Spring (February 28 - April 15) 54 

32 Linear correlation coefficients as an indicator of scatter in 

the data between emissivlty and API for Summer (April 16- June 8) . 55 

33 Linear regression slopes as an indicator of the sensitivity 

of the ESMR sensor to API for Fall (August 12-November 1) 56 

34 Linear regression slopes as an indicator of the sensitivity 

of the ESMR sensor to API for Winter (November 2-February 27). . . 57 

35 Linear regression slopes as an indicator of the sensitivity 

of the ESMR sensor to API for Spring (February 28-April 15). ... 58 

36 Linear regression slopes as an indicator of the sensitivity 

of the ESMR sensor to API for Summer (April 16-June 8) 59 

37 A plot of SSD (Model 1) and API for a low yield area 63 

38 A plot of SDD (Model 2) and API for a low yield area ....... 64 

39 A plot of SSD (Model 1) and API for a medium yield area 65 

iii 


LIST OF FIGURES (continued) 

Figure Page 

40 A plot of SSO (Model 2) and API for a medium yield area .66 

41 A plot of SOD (Model 1) and API for a high yield area 67 

42 A plot of SOO (Model 2) and API for a high yield area 68 

43 A time series plot of cumulative SOD (Model 1) for three 

selected grid points 69 

44 A time series plot of cumulative SOD (Model 2) for three 

selected grid points 70 

45 An example of the technique used to forecast yield 72 

46 A scatterplot of final SOI (Model 1) versus yield 73 

47 A scatterplot of final SOI (Model 2) versus yield. 74 

48 A graymap depicting the residuals from Model 1. The residuals 

were calculated as the difference between actual and predicted 
values of yield 76 

49 A graymap depicting the residuals from Model 2. The residuals 

were calculated as the difference between actual and predicted 
values of yield . 77 

50 A plot of cumulative SOD (Model 1) for three low yield grid 

points 78 

51 A plot of cumulative SOD (Model 1) for three medium yield 

grid points 79 

52 A plot of cumulative SOD (Model 1) for three high yield 

grid points 80 


LIST OF TABLES 


Table Page 

1 A comparison of CS weighting factors for Model 1. ....... . 61 

2 Comparison of SOD weighting schemes 81 


1 .0 INTRODUCTION 
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1.1 Necessity of This Research 

Estimates of surface soil moisture In a large area context are pr1> 
marlly useful for large area crop nenitoring, crop yield forecasting, 
estimating flood hazards^, and as Inputs Into (dynamic atmospheric 
models. Such estimates may also provide; Indications of soil moisture 
below the surface as well as provide a means for the determination of 
drought and the aerial extent of drought conditions. Conventional soil 
moisture measurements are very time consuming and not widely or regu- 
larly obtained over most of the United States. The spatial variations 
of soil moisture make It difficult to extrapolate conventional point 
measurements to represent an Integrated value over a large area. 

The benefits of crop yield forecasting are so wide ranging that 
some forecast methods are practiced In almost all countries (Idso ^ 
al ., 1975). A yield forecast Is of use In estimating supply, location 
of crop shortage areas, and In allocation of harvest resources such as 
harvesters and railroad cars. Inaccurate forecasts cause producers and 
shippers to make Inefficient use of facilities. Early prediction of 
crop yields would also allow adjustment of current Inventories. As more 
dependence Is placed on these forecasts, accuracy must be Improved. The 
loss In "social well-being" due to forecast errors has been modeled by 
Hayami and Peterson (1972). Their model computes the cost of a one per- 
cent error In a crop's yield forecast. When the model Is used with 1978 
United States data from Agricultural Statistics (1979) It Is calculated 
that a one percent error In wheat yield estimation would result In a 
$30,000,000 decrease In social well-being. This loss would be felt by 
consumers In the form of higher prices. 


The current crop yield forecasting technique used by the United 
States Department of Agriculture (USOA) includes the use of test plots. 
According to Idso (197S) ihe crop condition for this predictive 

estimate is visually determined. Cl imatologically predicted rainfall is 
then used as an additional input to forecast the final yield for the 
test plots. The results are extrapolated to the total acreage in the 
vicinity of the test plots. This method presents^ difficulties when 
applied on a large scale. Rainfall can be very localized, which may 
make the test plots unrepresentative of true local conditions. Wide 
variations in crop condition estimation can occur as the result of the 
subjective nature of visual observations. A large number of test 
observations must be used to obtain a statistically significant sample. 

The difficulty in using conmon meteorological variables and point 
measurements as inputs into soil moisture and crop models has led to 
research in remote sensing techniques to develop methods of obtaining 
Inexpensive and large scale assessments of soil moisture and crop condi- 
tion. Interpolation of point measurements and the resulting inaccura- 
cies would no longer be a probleei if the spatial variation can be 
mapped. Sensing devices, particularly when satellite mounted, give an 
integrated value for soil and crop conditions on a large scale; thus 
such a tool could be useful in several ways. Drought areas could be 
determined easily. During the growing season, the current crop condi- 
tions could be used to forecast a range of possible final yields. The 
final yield could also be predicted for inaccessible developing coun- 
tries and hostile areas. 
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1.2 Background 

The fundamental concepts of thermal microwave radiation that per- 
tain to this study are presented In detail In texts by Paris (1969 and 
1971). Newton (1977), Hess (1959), Marlon (1965), and Llntz and Simonett 
(1976). 


1.2.1 Soil Moisture Estimation 

The 1.55 cm ESMR Is a quasi -operational spacecraft system from 
which digital data can be acquired over the same area on approximately a 
3-day repeat cycle. This affords the opportunity to use time series 
data for multi -temporal mapping. Recent Investigations by Clhlar and 
Ulaby (1975), Meneely (1977), Schrougge et al . (1974, 1976a and 1976b), 
Schmugge (1976 and 1977) and Newton (1977) have demonstrated that sur- 
face emissivlty at the 1.55 cm. wavelength Is Inversely related to soil 
moisture ccntent In the surface layer. Sensitivity of this emissivlty 
to moisture content Is significantly diminished by an Increase In sur- 
face roughness, and/or an Increase In vegetation density. Consequently, 
the most significant results have been obtained on relatively bare, 
smooth soils. 

Schmugge et j^. (1977) presented case studies of ESMR*s spatial 
response to recent rainfall as related to vegetation and surface rough- 
ness. Relative vegetation densities were obtained from Landsat false 
color Infrared Images and surface roughness features were Inferred from 
U. S. Geological Survey surface land forms. Their results snow that 
variations In surface roughness and vegetative cover, plus the absence 
of large areas of bare soils, restrict the spatial mapping capabilities 
of soli moisture at satellite altitudes. 
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Temporal mapping of soil moisture shows a greater potential than 
other methods because the adverse effects of point-to-point variations 
of surface roughness and vegetation covir are minimized (McFarland and 
Blanchard, 1977). The moisture content and temperature of the emitting 
layer Integrated over the sensor footprint forms the major variations of 
temporal brightness temperature changes. Since the emitting layer 
temperature can be approximated, the brightness temperature changes may 
provide a fairly accurate Indication of soil moisture changes from rain- 
fall and subsequent drying. By using temporal mapping techniques with 
ESMR, McFarland and Blanchard obtained high correlations between micro- 
wave emissivlty and soil moisture modeled oy antecedent precipitation 
Indices (API) during the autumn (mlnimun vegetative period) over rela- 
tively flat terrain/ 

1.2.2 Plant Stress Estimation 

The Immediate effect of a water deficit Is a loss In plant turgor 
and closing of the stomates. This results In a rise In leaf temperature 
and a reduction In photosynthetic activity according to Ehrler and 
van Bavel (1967). When these effects occur the plant Is said to be 
undergoing stress. Chinoy (1962) studied the effects of a water deficit 
at several phenol ogical stages of the wheat plant. His results show 
stress acts to reduce the size and quantity of wheat grains. The reduc- 
tions were highest when the period of stress came In the latter stages 
of development. This Is due to the Increasing water requirement of the 
plant as It grows. A similar study by Day and Intalap (1970) found 
reduced grain yield and a lower weight per volume for wheat that had 
undergone stress. Day and Barmore (1971) showed that a stress period 
near heading lowered the quality of the resulting flour. 
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The use of highly Instrumented test plots has produced useful 
results. Ehrler (1973) measured plant canopy temperature (Ty) and 
ambient air temptrsture (T^). He found to be larger than Ty 
when the plant had adequate moisture supplies. As the soil dried out, 
the temperature difference became smaller. This study was expanded by 
Idso and Ehrler (1976). They found Ta to exceed Ty by up to 6*C for 
unstressed plants. For stressed plants, Ty exceeded Ta by up to 
8*C. The measurements were taken around local noon when the 
temperature differences should be near the maximum. 

A study of the diurnal change In leaf*a1r temperature differential 
was carried out by Ehrler et j^. (1978). They showed that Ty was 
larger than Ta for stressed plants during the daytime only. At night 
there was little difference In temperature differential between stressed 
and unstressed plants. They also showed that the diurnal change In 
temperature differential was twice as large for stressed plants as for 
unstressed plants. 

The difference between leaf temperature and ambient air temperature 
for stressed plants Is adaptable to remote sensing methods. Barthollc 
et al . (1972) used an airborne thermal scanner over cotton fields that 
were undergoing varying amounts of stress. Their results showed a high 
degree of correlation between leaf>a1r temperature differential and 
leaf-water potential. A study by Heilman ^ al^. (1976) developed a 
method of calculating evapotranspi ration using an airborne thermal 
scanner. The only drawback Is that the model requires ground level 
measurement of net radiation. 

The ability to use remote sensing techniques to obtain evapotran- 
spi ration values would be of use In models such as the one by Hller and 
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Clark (1971). They developed an Irrigation sc'.iedullng technique that 
Introduced the stress degree Index (SOI) defined as 

n 

SOI - ^ (SOi X CSi) (1) 

1-1 

The stress day factor (SO) Is characterized as: 

SO - I - (ETa/ETp) (2) 

where ETa Is the actual evapotranspi ration and ETp Is the potential 
evapotranspi ration. The crop susceptibility factor (CS) Is dependent 
upon the water requirement when stress Is encountered. It serves as a 
weighting function for the stress day factor. 

A study by Idso et a^. (1978) outlined the use of a stress degree 
day (SOD) summation. It Is analogous to the familiar growing degree day 
and Is formulated as 

n 

Yield « 3 - 4 SSOi (3) 

1-1 

where SOD Is the 2 pm value of 

SOO • Tygg - Tj^|» 

This model assumes a mixmum posslvle yield (a) with no stress 
encountered. The value of SOO Is changed to amount of yield reduction 
when multiplied by a. The model requires In situ measurement of air 
temperature (T^i^.) a few meters above the cc^op while the plant canopy 
temperature (Ty^g) can be remotely measured with an Infrared 


thermometer. 


1.3 Objectives 


The tMO satellite mounted microwave systems to be used were the 
Electrically Scanning Microwave Radiometer (ESMR) and the Scanning 
Multichannel Microwave Radiometer (St<1MR). This study only considered 
ESMR because processing difficulties of the SMMR data have delayed 
receipt of the data until early 1981. 

The thrust of this study was to use ESMR brightness temperatures 
(BT), meteorological data, climatological data and crop Information over 
the southern Great Plains to accomplish three objectives. 

The objectives are: 

1. Test and Improve correlations between passive microwave 
brightness temperature and Indices of soil moisture. 

2. Detect free water after storm frontal passage. 

3. Examine the feasibility of adapting a stress degree Index (SOI) 
that accepts passive microwave Input data In order to Indicate 
drought conditions and/or predict winter wheat yield. 

i 
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2.0 SITE DESCRIPTION 

A large region of the southern Great Plains was used In this study 
(Figure 1). This region was selected for two reasons: f1rs&, It encom- 

passes the area used by McFarland and Blanchard (1977) In the prelimin- 
ary study of ESMR data (Figure 2) » and, secondly, the area Is a prin- 
ciple source of hard winter wheat. The study region has a wide range of 
roughness, soil types, vegetative cover, and average yearly rainfall 
amounts. General characteristics cf the study region are: (1) there 

are no mountain ranges In this region; (2) the climate varies from arid 
to moderate; (3) the plains are cut by river headwater erosion; and (4) 
the eastern half of the preliminary study area has timber along the 
steams thickening along the eastern boundary to some solid timber 
areas. A general soils map for the preliminary study area Is shown In 
Figure 3. 

2.1 Data Available 

Basic data used In this study consisted of meteorological pol. 
measurements, passive microwave return of the ESMR, and crop yield 
Information. Meteorological data over the study area was obtained from 
the NOAA climatological network. The Information available Includes 
dally rainfall and maxlmun and mlnimun temperatures. These data lend 
well to gridding to 2S km because station spacing Is also approximately 
25 km. All data over the region was related to a 25 x 25 km grid. The 
block of grid points used by McFarland and shown In Figure 2 was used as 
a preliminary study area to test and Improve the relationship between 
passive microwave response and soil moisture Indices. 

ESMR brightness temperatures on digital magnetic tape were fur- 
nished by the Goddard Spaceflight Center for the period from August 1973 










ri<fure 3. Soils map of tne preliminary study area. 
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through May 1975. Data for the simmer months were not available. The 
Interval between days with ESMR data ranged between 1 and 7 days> with 
an average of approximately 3 days. The brightness temperatures used In 
the analysis were observed near 1700 GMT. 

Crop Information was derived from state agricultural statistics on 
a county basis. These data are very Inaccurate and become even more so 
when relaxed to the grid. Percentages of area devoted to winter wheat 
and total croplands by county are shown for the preliminary '..tudy area 
In Figure 4. The Inaccuracy of some of this Information becomes appar- 
ent when the special variability of winter wheat areas are observed from 
Landsat false color composites. 

The percentages of grid area covered by cropland and continuously 
cropped dryland winter wheat are shown In the form of graymaps, Figures 
5 and 6. Figure 7 Is a graymap of the 1974 dryland winter wheat yield 
(bushel s/ acre) . At the top of each graymap are the percentage limits 
for each gray shade. 

Several crops In addition to wheat are grown In the grid area. 
Since the coarse resolution of the ESMR precludes sensing only wheat 
fields, some contamination of the data Is Inescapable. However, a com- 
parison of Figures 5 and 6 shows that continuously cropped winter wheat 
Is the predominant crop In the central areas of the grid. 


2.2 Cropping Practices 

Normal cultural practices for winter wheat In Oklahoma and Texas 
begin the cropping season with plowing In July and seedbed preparation 
In late August to late September. Depending upon the availability of 
soil moisture, seeding takes place from mid-September through late 
October. The halfway point of planting Is near October 1. Early plant- 


















Figure 5. A graynap representing the percentage of each qridpoint 
covered by cropland. 














Fiqure 7. A graynap representing the 1974 continuously cropped dryland 
winter wheat yield, in bushels' per acre, for each grid point 
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ings to be used for winter pastures normally will have sufficient growth 
for grazing by December 1. Consequently, the soils are essentially bare 
from July into November for the winter wheat land of Oklahoma. 

Grazing begins around December 1 when the wheat Is 5 to 25 cm 
tall. Grazing continues until approximately two weeks before jointing, 
which occurs in March or April. The water used by the wheat crop 
Increases from about 0.23 cm per day In March to 0.89 cm per day In 
April. The wheat begins to boot In mid-April and continues to mature 
until harvest In June. 

Other crops that are grown In this area have different phenolo- 
gies. In the early Spring the crop> such as cotton and grain sorghums 
are being planted whereas the winter wheat Is already 15 to 25 cm tall. 
In May all crops and rangelands significant vegetative cover. The 
non-wheat crops continue to mature through the Summer. Most are harves- 
ted in October. 

The phenol ogical stage of the wheat crop varies up to four weeks 
due to latitudinal differences, with the southern crops leading in 
development. Figure 8 shows the number of weeks that wheat development 
in other areas lags that of the Texas crop. State agricultural statis- 
tics compiled by Whitehead (1979) were used as a basis to prepare 
Figure 8. 
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Figure 8. Nunber of weeks that the phenoiogicai development of wheat 
lags that of the Texas crop. 



3.0 NIMbJS 5 ESMR DATA 

The 5 satellite was launched In December 1972 Into a circu- 

lar orbit at an altitude of 1112 km. The satellite Is sun- synchronous 
with a local noun (ascending) and a midnight (descending) equator cros- 
sing. The period of the orbit Is about 107 min with successive orbits 
crossing the equator with longitudinal separation of 27 deg. This orbit 
affords daylight temporal coverage of the entire earth with approxi- 
mately a 3-day repeat cycle. One Instrument on Nimbus 5, the ESMR, Is a 
horizontally polarized radiometer with a central frequency of 19.35 GHz 
and an Intermediate frequency bandpass from 5 to 125 MHz; thus It Is 
sensitive to radiation from 19.22 to 19.475 GHz. Every 4 sec, the unit 
scans perpendicular to the spacecraft velocity vector from 50 deg to the 
left of nadir to 50 deg to the right of nadir In 78 discrete steps. 
From an altitude of 1112 km, the resolution Is 25 km x 25 Ion near nadir 
and degrades to 160 km crosstrack x 45 km downtrack at the ends of the 
scan, as shown In Figure 9 (Wllhelt, 1972). 

For this study antenna brightness temperatures (BT) were restricted 
to resolutions less than 50 km. Wllhelt (1973) reported a cross polar- 
ized grating lobe In the ESMR. This correction has not been Incorpor- 
ated In the data. 


175 H 
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Figure 9. ESMR ground resolution as a function of 
scan angle from nadir. 
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4.0 TECHNIQUES OF DATA HANDLING 

The brightness temperatures, precipitation, and air temperatures 
were resampled to emulate a grid having 25 km cells. The grid estab> 
11 shed Is based on a polar stereographic projection. 

The grid point values are determined from 

( 5 ) 

where c|in) Is the value of the variable at observation point m and 
Is the Interpolated value at grid point 1j. The distance-dependent 

weight function, w. Is a modified Barnes exponential (Barnes, 1973) 

* exp ( -a“) 

where d Is the distance from the grid point to the observation point and 
a Is the response parameter to control the shape of the weight function. 

4.1 The Soil Moisture Model 

In the absence of actual soil moisture measurements, a simple soil 
moisture model was used to account for changes In moisture In the ESMR 
emitting layer. An antecedent precipitation Index (API) was selected 
because of Its simplicity and Its ability to Infer upper-level soil 
moisture. The only Input required by the API model Is the precipitation 
which, for large areas. Is readily available at raingage locations from 
climatological data. Effective precipitation was considered a direct 
input to the soil water storage that Is estimated by the API. Losses of 
soil moisture due to evaporation and transpiration were assumed to 
decrease exponentially with time (Linsley, Kohler, and Paulhus, 1975). 
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Shown mathematically, the relationship for a single rainfall event Is: 

APl<+t • ( 7 ) 

where P Is effective precipitation, 1 Is the day number when rainfall 
occurred, t Is the time after rainfall, and k is a recession factor 
which accounts 'or seasonal differences In evapo transpiration losses. 

Rather than total the combined Influence of all the rainfall events 
In a period, dally Indexes can be calculated by setting t equal to 1 In 
equation 7 and repeatedly using the following form (Saxton and Lenz, 
1967): 

APIi » Pi + (APIi.i 

Before the first API value was used In a correlation with emlssl- 
vlty, the API model was allowed to stabilize by using a mlnlmun of 30 
days of rainfall history. The relationship between rainfall amount (R) 
and effective precipitation, developed by Blanchard et (1979), was 
used to account for runoff. Am empirical recession curve developed at 
the SEA/AR (Southern Great Plains Watershed Research Center, Chlckasha, 
Oklahoma) by OeCoursey (1974) was used for calculating the dally reces- 
sion factor, k. The final form of the soil moisture model was 

891 

APIi * Rl (APIi-i * ki) (9) 

where 

Ri * dally rainfall amount (cm). 

4.2 Emissivlty Model 

The temperature of the emitting layer was approximated by the dally 
maximum temperature (Tmt) emissivlty model, 

c * Tbt/Tmt (10) 

where Tbt Is the ESMR brightness temperature and c Is the emissivlty. 


The emitting layer for the short wavelength ESMR Is most likely limited 
to the top centimeter of the surface being observed. The overpass time 
of ESMR over the study area Is near local noon and maximum air tempera- 
ture usually occurs several hours later. Because the maxlmun soil 
temperature usually leads the maximum air temperature, this Is believed 
to be a sound approximation. The sensitivity of this model to errors In 
the emitting layer temperature Is small. An error of 10*K will only 
yield a 4 percent error In the predicted emissivlty. 

Since the study period encompassed winter periods, consideration 
was made of the Influence of frozen ground on the ESMR return. The 
dielectric properties of Ice are completely different from water. Ice 
Is not a dipole molecule and has a low dielectric constant. When a soil 
Is frozen the emissivlty Is high and Independent of soil moisture. 
Therefore, to avoid the confusion possible from observations of frozen 
ground, the emissivlty values were not used In this study If the maximum 
air temperature was less than 283**K. 

Rain in an atmosphere can also have a very significant effect on 
upwelling radiation, Paris (1971). Radar summary charts from the NWS 
were used In determining the presence of rain between the ESMR and the 
earth surface for each grid point. When rainfall occurred at the time 
of the ESMR overpass, the emissivlty data for affected grid points were 
omitted. 

4.3 Stress Degree Index Model 

This study exploited the difference between canopy temperature and 
air temperature for stress measurement. Two models were used to 
describe stress. One model requires only remotely sensed data and 
climatological temperatures. The other model uses air temperatures 


modeled from dally maximum and minimum temperatures. These two models 
were chosen so that this stu<^ could not only determine whether or not 
the stress measurement works, but also whether or not local temperature 
data were needed to make It work. 

The general form of the stress degree Index (SOI) model used In 
this study was patterned after the Hller and Clark (1971) model. The 
form of the model used was 

n 

SOI - I (SOOi X CSi) 

1«1 

where SOOi * number of stress degree days for day 1, 
and CSi * susceptibility factor for day 1, 
n > number of days In the growing season. 

Both models calculated SOI using 

SOOi « tveg - Tair If Tveg^Tair 

or (12) 

SOOi , 0 If Tveg<TaIR 

Tveg f*’® vegetation temperature and Tair f^® temperature 
at the time of overpass. Tveg calculated using 

Tveg * V0*92 (i3) 

The 0.92 was assumed to be an average emissivlty of a vegetated 
surface. This number was chosen because It Is a representative value of 
emissivlty when the scene Is dry and/or vegetated. The first model, 
hereafter called Mod<tl 1, defined Tair as 

Tair “ 0.75 (tmax - Tmin) + Tmin. 

where Tmax ^)^® <la^ly maximum temperature and Tmin is the dally 
minimum temperature. Tair estimation of the air temperature at 

the time of the 11 a.m. ESMR sensor overpass. 
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Both models calculated SOI using 


SDOi * Tveg - Tair If TvEG^UlR 

or (12) 

SOOi . 0 If Tv£6<T*m 

Tveg vegetation temperature and Tair Is the air temperature 

at the time of overpass. Tveg calculated using 

Tveg “ Tb/o.92 (13) 

The 0.92 was assumed to be an average emissivlty of a vegetated 
surface. This number was chosen because It Is a representative value of 
emissivlty when the scene Is dry and/or vegetated. The first model, 
hereafter called Model 1, defined Tair as 

Tair « 0.75 (Tmax - Tmin) + Thin. (h) 

where T^ax ^s the dally maximum temperature and T|ti|xi«| Is the dally 
minimum temperature. Tair ^s an estimation of the air temperature at 
the time of the 11 a.m. ESMR sensor overpass. 

The CS curve used with Model 1 Is shown In Figure 10. The stan- 
dardized CS curve was calculated using water use requirements measured 
for winter wheat at the Bushland, Texas USOA research station. Since 
phenological development Is slower In the northern portion of the grid, 
the CS curve must shift to the right when used north of the Texas area. 
The amount of the shift. In weeks. Is shown for each grid point In 
Figure 8. 
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•utcspUbiUiy curve developed lor tide 


The second model, hereafter referred to as Model 2. required only 
microwave brightness temperature as a real time Input. The SOD was for* 
mulated as before only Tavg climatological estimate of the 

temperature at the time of overpass. It was calculated using the 
average dally maximum and minimum temperatures for Oklahoma City, 
Oklahoma and Wichita, Kansas. These two cities were chosen because they 
He In a predominantly winter wheat area. For each day of the summation 
the average maximum temperatures for the two cities were averaged. The 
same was done for each day's minimum temperature. The resulting two 
curves are shown In Figure 11. From these curves a climatological 
estimate of maximum and mlnlmun air temperature for the central grid 
area can be obtained. These estimates were then used with Eq. (14) to 
calculated Tair. The remaining computation proceeded as In Model 1. 
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5.0 AN/>4.YSIS AND DISCUSSION 


5 . 1 Correlation Between API and Em1ss1v1t>> 

This section was directed toward developing relationships between 
einissivlty and API over the training grid area In Oklahoma (Figure 2 ), 
The results were presented In a thesis by The Is (1979) and In a paper 
presented at the 1979 AGU Spring Meeting (abstracts In Appendix A). 
The results are summarized and presented here. 

The calendar year was divided Into four new-standard ESMft seasons 
because of climatological factors, crop phe>io1og1es and cultivation 
practices. Percentages of areas devoted to winter wheat and total crop- 
lands for the grid area are presented In Figure 4. Fall was defined 
from August 12 to November 1 which corresponded to a mlnlmun vegetation 
i:er1od when fields are relatively flat. This Is the period studied by 
McFarland and Blanchard (1977) and their results were duplicated. The 
fall correlation coefficients between emissivlty and API for each grid 
point are presented In Figure 12. By comparing the areas of greater 
than 0.60 to Figure 4 and Landsat color composites (not presented) It Is 
apparent that cultivated agricultural lands give the best correlations. 

Winter, as defined by this study (November 2-February 27), Is char- 
acterized by periods of frozen soil surfaces and upward movement of 
moisture due to temperature gradients. The correlations during this 
period are much less significant with values generally around -0.6C. 
Both the emissivlty and API models are not well suited for the winter. 
The API model Is very simple and does not account for movement of water 
due to temperature gradients. The dielectric properties of Ice are 
significantly different from water, making the emissivlty Independent of 
soil moisture In the frozen state. 
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Spring (February 28>Apr11 15) was defined as a rather short season 
generally bounded by the end of frozen soil surfaces and the beginning 
of the winter wheat "boot" stage. This Is a period of smooth soil 
surface with Increasing vegetative cover over winter wheat acreage. 
Over other agricultural croplands the soil surface Is bare, but Is 
either In a rough bedded condition or freshly planted In rows. The good 
correlation areas, shown In Figure 13, generally correspond with 
predominantly winter wheat agricultural areas. This shows up quite well 
when the grid Is overlayed onto a Landsat color composite taken during 
April when the only growing vegetation Is winter wheat. 

All croplands are densely vegetated during summer (April 16-June 8) 
as defined In this study. Correlation coefficients are corresponding 
poorly with values averaging around -0.50. 

The differences between the spring correlations of winter wheat and 
non-wheat croplands were Investigated further by plotting six grid 
points from each area. The plots for fall, shown In Figures 14 and 15, 
Indicate very little differences In the two areas. The values of the 
slopes and Intercept agree closely with those obtained by McFarland and 
Blanchard (1977) (slope » -0.0232, Intercept » 0.92). It should be 
pointed out that McFarland used one year's data and this study used two. 

Spring scatter plots are presented In Figure 16 and 17. Correla- 
tions for the winter wheat area are significantly higher than non-wheat 
croplands. The slope for the winter wheat areas has slightly 

decreased from the fall value (-0.0227 to -0.0156). This Indicates that 
the small winter wheat vegetation may affect but does not destroy the 
good relationship. The suniner slope for the same grid points decreases 
to -0.0081 with a corresponding correlation coefficient of -0.48. Our- 
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NON-WHEflT CROPLANDS FALL (AUG12-N0V1I 



CORR. COFFF. = -0.8333 * SLOPE = -0.0213 





HINTER WHEAT AREAS SPRING (FEB28-APR1 51 






CORR. COEFF. = -0.5835 " SLOPE = -0.0093 

Figure 17. Spring scatter plot for six non-wheat croplands grid points. 




Ing summer the vegetation has reached a threshold density so that the 
ESMR's response to API becomes masked by the vegetation. 

5.2 Free Water After Storm Frontal Passage 

The same techniques were used to map both API and emissivlty over a 
large portion of the Southern Great Plains. 

API and emissivlty gray maps were made for each day that an ESMR 
pass occurred over the test region. API values were calculated for each 
day In each cell and then merged Into a spatial map for each day. 
Density slices were chosen for the API gray scale such that Intermediate 
steps In gray tone represent a range of one centimeter while the maximum 
and minimum gray tone may represent more or less than one centimeter. 
Gray tones were selected for emissivlty maps by first fitting a regres- 
sion line to the API and emissivlty data for each date and calculating 
emissivlty values to correspond to the API divisions. This admittedly 
enhances the visual correlation between the two gray scale maps for each 
date. It Is significant, however, that differences In the scales for 
emissivlty are minor throughout the different seasons. 

Figures 18 and 19, 20 and 21, 22 and 23, 24 and 25, 26 and 27 are 
pairs of gray scale maps that have been selected to Illustrate moisture 
conditions after major storms. In each pair of maps for a particular 
date, a remarkable correlation can be seen between the ESMR emissivlty 
and the API. 

A preliminary 16mm film was produced with dally gray maps of API 
and emissivlty on each frame. By visual Inspection, several model Ifig 
problems appeared. Whenever the dynamic range of API was small (low API 
throughout the grid) correlation was poor. When the range of API was 
large there was a good correlation. API grey maps during the winter 
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EMISSI VITT 


DRTE: 3/13/75 


MERN = 0. 878 STRNORRD OE V I RT I 0N = 0 . 042 
0.772 0.833 0.845 0.856 0.868 0.880 0.891 0.903 0.955 






months indicate that there is insufficient depletion in the API model 
after a significant rainfall event. The winter time recession factor 
(k) was developed fo; soil moisture to a depth of 9 in. and it appears 
that values of k are too high during the winter period. The values of k 
may also have a latitudinal dependence when applied to the soil surface. 

Correlation coefficients for each cell were calculated after divi- 
ding the entire data set Into four parts to represent major portions of 
the growth cycle for winter wheat. The planting, emergence and stool Ing 
of the plant take place In the fall. August 12 to November 1 was 
selected as an appropriate time period for what we define as the fall 
period. Vegetation during this period seldom exceeds 15 centimeters In 
height and large percentages of bare soil are normal. Winter, as 
defined for this study, extends from November 2 through February 27. 
The winter season Is characterized by periods of frozen ground when 
upward movement of moisture due to temperature gradients can be expected 
and downward movement or Infiltration may be restricted when rain 
occurs. 

The spring period from February 28 to April 15 is characterized by 
rapid growth to the booting stage of wheat. In areas where heavy winter 
grazing Is practical In Oklahoma and Texas, It Is common to remove 
cattle prior to April 15 to prevent loss In the grain crop. Near the 
northern edge of Kansas, the boot stage for winter wheat may be two to 
three weeks later, while this stage may be reached In the south end of 
the study area as early as April 1. After April 15 to the time of wheat 
harvest has been defined as early sunmer. During all of this period, 
the grain crop is forming and a minimum of bare soil Is visible until 
the plant begins senescence. 
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A map showing the distribution of wheat crop as a percentage of each 
cell, Figure 28, provides an overview of the agricultural distribution 
within the study area. This map points out the portions of the area 
that are of significant economic Importance If we wish to predict 
moisture availability for wheat yield estimates. Figures 29, 30, 31 and 
32 Illustrate the fact that the portions of the study devoted to wheat 
show better correlation between the emissivlty and API values during the 
fall season. Figure 30. An examination of all four figures depicting 
correlation Indicates that there Is always low correlation In the 
eastern central portion of Oklahoma. This area of poor correlation Is 
recognised as the moderately rough timbered area typical of the cross 
timber region. In general, these Illustrations Indicate that the rela- 
tively flat farming areas of the great plains lend themselves to the use 
of the ESMR system as a moisture monitoring tool. 

Another measure of the usefulness of the ESMR system can be deter- 
mined from examination of the slope of a regression Tine relating 
emissivlty to the API value within each cell. Grouping the data In the 
sanie four periods used before, regression slopes were calculated for 
each cell and the slope was tnen scaled to produce the grayscale maps 
shown In Figures 33, 34, 35 and 36. These figures show that sensitivity 
of the system to changes In API tends to be greatest In the northwest 
quadrant of the study area. Again, the slope appears to be even more 
closely related to the smoothness of the terrain. The lowest sensitiv- 
ity, as well as the poorest correlations, was found In the winter 
period. Fortunately, this period seldom has drought severe enough to 
endanger wheat yield. 
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Figure 20. Percent of each grid area covered by winter wheat (derived by 
state agricultural statistics). 
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Figure 29. Linear correlation coefficients as an indicator of 
scatter in the data between emissivity and API for 
Fall (August, 12-November 1). 
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Linear correlation coefficients as an indicator of scatter 
in the data between emissivity and API for Sprinq 
(February 28-April 15). 
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Figure 32. Linear correlation coefficients as an indicator of scatter 
in the data between emissivity and API for Summer (April 16- 
June 8) . 
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Figure 33. Linear regression slopes as an indicator of the sensitivity 
of the ESMR sensor to API for Fall (August 12-‘'ovember 1).’ 
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Linear regression slopes as an indicator of ,.he sensitivity 
of the ESMR sensor to API for Winter (November 2- 
February 27). 
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Figure 35. Linear regression slopes as an indicator of the sensitivity 
of the ESMR sensor to API for Spring (February 28-Apri1 15). 
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Figure 36. Linear regression slopes as an indicator of the sensitivity 
of the ESMR sensor to API for Summer (April 16-June 8). 
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5.3 Drought Condit'lons and Winter Wheat Yield 


The results of this section were presented In a thesis by Richter 
(1980) (abstract In Appendix A). 

As this study of stress Indications from the ESMR data progressed, 
some restrictions were made that limited the area used for final 
analysis. Grid points that did not represent areas of at least 30 
percent winter wheat were not used. This was to Insure that wheat 
emission made a significant contribution to the brightness temperature. 
A further restriction was Imposed that the wheat fields comprising the 
30 percent must be continuously cropped. This removed from 
consideration areas of simmer fallowed and Irrigated wheat land. Both 
of these cropping methods Increase plant available moisture In 
comparison to continuous cropping that depends on rainfall In the same 
year the crop Is grown. A total of 95 grid points fit the above 
criteria. 

Contrary to expectations. It was found for both models that the 
crop susceptibility curve acted to decrease model accuracy. The compu- 
ted curve In Figure 10 has a factor of 10 difference between maxlmim and 
minimum values. The shape of the curve was vertically compressed to 
study the effect of a smaller range In susceptibility. This was 
accomplished by talcing the CS to an exponent before multiplication with 
the SOD. The exponents that were tested are listed In Table 1 along 
with the resulting correlation coefficients. The highest overall 
correlations between SOI and winter wheat yield resulted when the curve 
was compressed Into a horizontal line. In effect, this study found that 
all springtime stress had an equal effect on yield reduction. This is 
In contradiction to the findings that a stress period near heading 
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TABLE 1 


where 


A comparison of CS weighting factors for Model 1 


QSk 

Correlation Coefficient 

X • 1.0 (F1g.6) 

- 0.17 

• 1.5 

- 0.10 

■ 0.4 

- 0.35 

• 0.2 

- 0.43 

- 0.03 

- 0.49 

■ 0.0 

- 0.51 


reduces yield more than If the period came earlier (Chlnoy 1962). 
HoMever, other studies of Mheat by Robins and Domingo (1962) have found 
little difference In CS magnitude between different phenol ogical 
stages. The Importance of early stress Is related to the storage of 

soil moisture. If a cron receives abundant moisture In the early 
months. It not only reduces stress then but also throughout the season. 

Since the CS curve acted to reduce the Importance of this period, the 

curve also r<c-Juced the accuracy of the yield estimate. Because of this, 
a value of I.O was used for CS at all times. This has the effect of not 
using a CS curve at all. Figures 37-42 are time series plots of SOI and 
API grid points that differ In yield. The Inverse relationship between 
SOI and API Is quite evident. A close examination shows that If the API 
Is greater than 3.0 cm, then neither model shows stress occurring. 

After a large rainfall event, with an API above 8 cm, the models do not 
show stress until the API decays to around 2.0 cm. In this study dally 
values of TygQ were obtained by linear Interpolation of Tyge between 
the days of ESMR data. This leads to a dependence of SOI on air 

temperature because the Tv£g values could not fluctuate with air 

temperature when the values were Interpolated. 

A seasonal summation of weekly SOD totals for the models Is shown 

In Figs. 43 and 44. A representative gridpoint from each of three 

categories, I.e., low, medlun and high yield, was chosen for the plots. 
The figures show additional justification for not using the CS curve In 
Fig. 10. Note that much of the yield reducing stress came In the first 
three months. The Importance of this period Is greatly decreased when 
the CS curve Is used. 
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SDI MODEL 1 
GRID POINT 11117 
SOI = V29 

TIELD = m tBU./fiC.) 
SMOOTH PLAINS 












SDI MODEL t 
GRID POINT SU 
SDI = 377 

riELO = 20 (BU./flC.) 
IRREGULAR PLAINS 













Figure 41. A plot of SOD (Model 1) 


and API for a high yield 


SOI MODEL 2 
GRID POINT 716 
SOI - 286 

YIELD = 32 (BU./BC.) 
IRREGULAR PLAINS 







Figure 43. A time series plot of cumulative SDD (Model 1) for three selected grid pot 






Figs. 43 and 44 are good illustrations of the yield prediction 
capability of the models. The yields are forecast In the following 
manner. First, within areas of homogeneous roughness, all SOI totals 
for a particular yield are averaged. The yields would then be plotted 
on the right-hand side of the figure across from the corresponding 
average SOI. For this study, gridpoints In areas of Irregular plains 
were used to compute the averages. The result Is shown In Fig. 45. The 
figure also shows an example of yield prediction. The current date and 
SOD total are located on the graph. The curves from Fig. 43 
representing maxlmim and mlnlmun stress, with yields of 10 and 30 bu/ac 
respectively, are vertically shifted until they pass through the point. 
The Intersection of these curves with the yield axis gives an estimate 
of maximum and minimum yields. From the example In Fig. 45, the 
estimate would be between 18 and 22 bu/ac. 

The task of predicting most probable yield is much more difficult. 
While It would certainly be between the maxlmun and minimum forecast 
yields, the exact location Is difficult to predict. This problem could 
best be solved by the Incorporation of a meteorological forecast. For 
example. If the long range forecast predicted above average temperatures 
and below normal precipitation, an above average amount of plant stress 
would be expected. If this were the case, the probable yield would be 
closer to the maximum stress estimate of yield rather than the minimum 
stress estimate. Thus, the forecast of probable yield would then be 
adjusted with the amount of adjustment dependent upon the forecasted 
departure from normal weather conditions. 

Figs. 46 and 47 are plots showing the variation between SOI and 
yield for the 95 grid points. Correlation coefficients, as well as the 
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Figure 46. A scatterplot of final SDI (model 1) versus yield. 




Figure 








slope and Intercept, are shown In each figure. It is interesting to 
note that Model Z has a higher correlation coefficient than Model 1. 
This was unexpected because Model 2 uses climatological air temperatures 
instead of the measured temperatures used in Model 1. The encouraging 
results obtained with Model 1. indicate that a yield forecasting 
technique using only remotely sensed data is feasible. 

Figs. 48 and 49 are graymaps showing the residuals calculated by 
linear regression. lT»e slopes and intercepts shown in Figs. 46 and 47 
were used to calculate the residuals. For each grid, residuals were 
compared to the percentage of continuously cropped winter wheat. There 
was no decernable pattern when residuals were plotted against percent 
winter wheat. The accuracy of the models was therefore unresponsive to 
differences in the grid point's wheat content. 

The effect of surface roughness is shown in Figs. 50 - 52. For 

each of three yields, three gridpoints were chosen that represented the 
widest range of final SOI. The gridpoints were then located on a map of 
uses classification of land form. This map outlined areas of smooth 
plains, irregular plains, tableland, low hills, and mountains. As the 
figures show, the magnitude of the SOI is directly related to the 

roughness of the land. Thus, much of the scatter in Figs. 46 and 47 is 
due to variations in roughness. 

An investigation was undertaken to examine the method of SOD 

summation. That is, do ten days of one SDO per day have the same yield 
reducing effect as one day of ten SOOs? To answer this question, three 
additional SCO summations were made. Each daily SDO value was taken to 
exponents of 1.3, 0.7, and 2.0. Table 2 shows the effect of the 

exponents on the summation and the resulting effect on the correlation 
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Figure 48. A graymap depicting the residuals from Model 1. The 
residuals were calculated as the difference between 
actual and predicted values of yield. 
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Figure 49. A graymap depicting the residuals from Model 2. The 
residuals were calculated as the difference between 
actual and predicted values of yield. 
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6.0 CONCLUSION 

Agricultural lands show the greatest potential to use the 
Electrically Scanning Microwave Radiometer to Infer soil moisture. 
These usually are better soils situated on smoother, less hilly or less 
eroded land. Tillage practices are such as to afford periods of smooth 
and bare soils during the year. In contrast, pastures, rangelands and 
cross timber region usually are situated on poorer soils (some rocky) 
and the terrain which Is unsuitable for agriculture. The surface of 

untllled lands Is almost never completely bare. As a rule. It Is 
covered with growing or dead vegetation or with timber. 

Results over the predominant winter wheat areas Indicate that the 
best potential to Infer soil moisture occurs during fall and spring. 
These periods encompass the growth stages when soil moisture Is most 
Important to winter wheat's yield. 

This study has also shown that It Is possible to Identify wheat 
stress periods with the ESMR using the stress degree day approach. 
However, the amount of scatter Indicates that the models need further 
development. The scatter In the results Is attributable to several 
sources. The effect of variations In the look angle Is noticeable In the 
data. There were Instances when the ESMR sensed the same grid point on 
consecutive days. Usually one of the brightness temperatures was much 
too low. Upon close examination. It was found that the anomalously low 
temperature was usually sensed at the extremes of the scan angle. The 
major sources of Inaccuracy are the subjective yield estimates used for 
comparison with the SOI. These estimates gave one value of yield for an 
entire county, which often contained several grid points. Thus, one 
Inaccurate county estimate can cause high residuals for several grid 


points. A possible example of this Is the two grid points with a 10 
bu./ac. yield estimate. Both points used the same county for a yield 
estimate. The models predicted yields of 16 to 20 bu./ac., similar to 
yields In surrounding counties. These two points had the highest 
residuals for both models. 

Due to the avallablity of only one year of data, analysis of stress 
detection was limited. A multi-year study Is needed to validate these 
results. This would allow the calibration of each grid point 

separately, minimizing point to point variations. Graphs such as those 
In Figs. 43 and 44 could be compiled for each grid point. This should 
allow for more accurate forecasts. 

The model which used climatological air temperatures was more 

accurate than the model that used actual temperatures. This Is probably 

due to a chance variation and cannot be expected to hold true for 

succeeding years. However, the high correlation coefficient obtained 
with Model 2 Indicated that a good forecast can be made using only 
remotely sensed data. Perhaps the best combination would be 

simultaneous measurement by a microwave radiometer and a thermal 
scanner. In this way. Inaccuracies due to the air temperature 

estimation can be eliminated. 

While this study concentrated on winter wheat, the SOI concept 
could be applied to other crops. However, the crop must meet certain 
criteria. The vegetation must be dense, so the sensor measures only 
plant canopy temperature. Sorghum, corn and rangeland grasses are 
likely candidates for this type of monitoring. For reasons explained 
earlier, two further restrictions must be remembered. First, the crop 
cannot have a moisture supply that Is unavailable to the surrounding 
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area. Second, the crop must fill a significant proportion of the grid 
cell. 

The results of this research clearly demonstrate the feasibility of 
using the short wavelength ESMR (1) to Infer soil moisture when surface 
conditions are acceptable. (2) to map the spatial distribution of 
moisture from storms resulting from frontal passages, and, (3) with 
further refinement, to develop an alorlthn to Identify drought 
conditions and to forecast crop yield. Atmospheric contributions, 
vegetation, and roughness produce detrimental effects to soil moisture 
detection at short microwave wavelengths. Longer wavelength radiometers 
such as the Scanning Mul t1 frequence Microwave Radiometer (SMMR) should 
lessen these effects when they become available. 

The ESMR or a similar system might be used operationally In the 
Interim period before longer wavelength passive microwave sensors are 
available. 
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APPENDIX A 


ABSTRACT 

CORRELATION OF BRIGHTNESS TEMPERATURES FROM THE ELECTRICALLY SCANNING 
MICROWAVE RADIOMETER (ESMR) WITH ANTECEDENT PRECIPITATION INDICES (API) 

M. J. McFarland (Dept, of Agricultural Engineering, 

Texas ASM University, College Station, Tx.) 

B. J. Blanchard, S. W. Thels (Remote Sensing Center, 

Texas ASM University, College Station, Tx.) 

Estimates of soil moisture In a large area context are primarily 
useful for large area crop monitoring and for estimation flood hazards 
on large and small drainage areas. Such estimates may also provide 
Indications of soil moisture below the surface as well as provide a 

means for the determination of drought and areal extent of drought 

conditions. A preliminary study correlated digital data from the ESMR 

to API that In turn Is correlated to soil moisture over the northwestern 

third of Oklahoma. Encouraging results were obtained for a three month 
period In the fall of 1973 when vegetation was sparse. Since the ESMR 
Is a quasi ’•operational system, data can be acquired over the same area 
on a three day repeat cycle. This provides the opportunity to 
Investigate changes In soil moisture through a time series. Temporal 
mapping reduces the effects of the point to point variations In 
vegetative cover, surface roughness, and soil characteristics. 

The moisture content and temperature of the emitting layer 
Integrated over the sensor footprint forms the major variations of the 
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temporal brightness temperature changes. Since the emitting layer 
temperature can be approximated, the brightness temperature changes may 
provide a fairly accurate Indication of soil moisture changes for 
rainfall and subsequent drying. 


ABSTRACT 


ESTIMATION OF SOIL MOISTURE WITH API ALGORITHMS 
AND MICROWAVE EMISSION 

Bruce J. Blanchard, Marshall J. McFarland 
Thomas J. Schmugge, and Edd Rhoades^ 

Large area soil moisture estimations are required for global sys- 
tems of crop yield estimation and flood prediction. Microwave sensor 

systems that as yet can only detect moisture at the surface have been 
suggested as a means of acquiring large area estimates. Measurements of 
soil moisture were studied to understand the correlation and Intercorre- 
latlon between moisture In surface soil layers and moisture In deeper 
layers. Relations previously discovered between microwave emission at 
the 1.55 cm. wavelength and surface moisture as represented by an ante- 
cedent precipitation Index were used to provide a pseudo Infiltration 
estimation. Infiltration estimation based on surface wetness estimated 
on a dally basis were used to estimate soil moisture at a depth of 15 

cm. by use of a modified antecedent precipitation Index with good 

results (r 2 * .‘/OlO and r 2 = .7383). The technique was modified and 


^Respectively, Remote Sensing Center, Texas A4M University, 
College Station, Texas 77843; Remote Sensing Center, Texas A4M 
University, College Station, Texas 77843, NASA/Goddard Laboratory for 
Atmospheric Sciences, Goddard Space Flight Center, Greenbelt, Maryland 
20771; SEA-AR, Chickasha, Oklahoma. 
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used to estimate soil moisture at 15 cm. depth when only an estimate of 
surface moisture each three days was available. Predictions based on 
estimation of surface wetness at three day Intervals resulted In r 2 
value of .6811 and J076 for the same date sets. The algorithms 
developed In this study can be used over relatively flat agricultural 
Liinds to provide Improved estimates of soil moisture to a depth greater 
than the depth of penetration for the sensor. 
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ABSTRACT 


Wheat Stress Measurement with the 
Electrically Scanning Microwave Radiometer (ESMR). (December 1980) 
John Charles Richter, B. S., Texas A4M University 
Chairman of Advisory Committee: Or. George Huebner 

This study has examined the feasibility of using a stress degree 
Index (SDI) to predict winter wheat yield. Microwave brightness 
temperatures In the spring of 1974 were used as an Input. Two 
formulations of the SDI model were studied. The difference was In the 
amount of real time data needed for the computation. The first method 
used only microwave brightness temperatures as sensed by the N1mbus>5 
ESMR. The second method also required data to estimate air temperature 
at the time of satellite overpass. 

Good correlations were found between winter wheat yield and SDI for 
both models. Forecast Inaccuracies were attributed to several sources. 
The primary source Is roughness variations between gridpoints. Other 
sources Included Inaccurate county agricultural statistics and 
differences In wheat content among gridpoints. 


